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Abstract. We study a Helmholtz-type spectral problem related to the prop- 
agation of electromagnetic waves in photonic crystal waveguides. The waveg- 
uide is created by introducing a linear defect into a two-dimensional periodic 
medium. The defect is infinitely extended and aligned with one of the coordi- 
nate axes. The perturbation is expected to introduce guided mode spectrum 
inside the band gaps of the fully periodic, unperturbed spectral problem. In the 
first part of the paper, we prove that, somewhat unexpectedly, guided mode 
spectrum can be created by arbitrarily "small" perturbations. Secondly we 
show that, after performing a Floquet decomposition in the axial direction of 
the waveguide, for any fixed value of the quasi-momentum k x the perturbation 
generates at most finitely many new eigenvalues inside the gap. 



1. Introduction 

The concept of a photonic crystal (also called a photonic band-gap material) was 
suggested in 1987 (see e.g. [9] for a textbook introduction) and has received signif- 
icant attention from both the theoretical and experimental viewpoint. In practice, 
photonic crystals are often manufactured using periodic crystalline structures, a 
feature of which is their ability to allow, or deny the propagation of electromag- 
netic radiation which lies in a well defined range of the frequency spectrum. We 
call ranges where the electromagnetic radiation can propagate spectral bands and 
ranges where propagation is prevented spectral gaps. These concepts will be made 
precise in Section [2] 

One possible application of photonic crystals is their use for manufacturing highly 
efficient optical waveguides which allow propagation of electromagnetic radiation 
only of very focussed frequencies. Briefly, these are created by taking some photonic 
band-gap material, called the bulk, and introducing a linear defect which breaches 
the periodicity. This may have the effect of allowing propagation of electromagnetic 
radiation in a range of the frequency spectrum in which propagation is not possible 
in the bulk. Waves of these frequencies, sometimes called guided modes, should 
then be highly focussed and almost entirely confined to the defect. 

An appropriate mathematical model for such materials is given by the spectral 
problem for the Maxwell equations, or in the context of polarised waves in two 
dimensions, by the Helmholtz equation. In this paper, we shall study the spectral 
problem for 

in R 2 , where e is the dielectric function (or, equivalently, the square of the refractive 
index of the material) . We wish to model a waveguide in a periodic crystal and to 
do this we take 

£ = £o + £l 
1 
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where eq is periodic and ei is a perturbation supported only in the waveguide, 
which we choose as a strip in the x-direction. 

In order to establish the existence of guided modes in a band gap of the un- 
perturbed medium (modelled by —Sq 1 A), we must prove that this gap of the un- 
perturbed problem now contains spectrum of — e _1 A induced by the perturbation 
e±. The mathematical investigation of these matters was begun in p] [21 Q3J [15] 
and [M], the latter two articles working with the full 3D Maxwell equations. In 
particular, [T21 HH [TS] give sufficient conditions for the existence of guided modes 
in the spectral gap of the unperturbed problem, for example: sufficiently wide gaps 
always contain spectrum of the perturbed problem. We note that existence of gaps 
in the spectrum of some problems with periodic background media was proved in 
PJ [3] and in [8] for the full Maxwell case. 

In this paper, we take another point of view and ask the following question: first 
fix a certain gap, no matter how wide, is there a certain threshold strength of the 
perturbation needed to produce spectrum in the gap? It turns out that spectrum 
appears in the gap of the unperturbed problem under arbitrarily small perturba- 
tions. We then further proceed to study the structure of the induced spectrum 
for fixed quasimomentum k x in the direction of the guided wave propagation. In 
particular, we show that for fixed k x , the eigenvalues introduced into the gap by 
the perturbation do not accumulate at the ends of the gap interval. 

It is interesting to compare our results to those available in the literature. In 
[3 [6] the authors show generation of a finite number of eigenvalues in the gaps 
by localized perturbations of the periodic medium. Another of their results in 
dimension d — 3 states that perturbations with E\ < (positive defects in their 
terminology) need a certain threshold strength to create eigenvalues in gaps. This 
is in sharp contrast to our result, which uses specific features of the dimension two 
nature of the problem; moreover, our defects are not localized. In [H Corollary 1], 
a result concerning the existence of a threshold for a two-dimensional waveguide 
analogous to that in [6] is claimed. This seems to be at odds with the result of 
Theorem [T] presented here. 

Our paper is structured as follows: in section[2l we introduce the operators to be 
studied and remind the reader briefly of the Floquet transform. Section [3] contains 
some preparatory material on the band functions and Bloch functions, which will 
play an important role later on. In section [4j we use variational arguments to 
prove the existence of guided mode spectrum. Finally, in section [SJ we consider the 
question of non-accumulation of eigenvalues at the ends of the gap. 



2. Formulation of the problem 

We first consider the fully periodic background problem in M. 2 . Let Eq be a 
positive bounded function on M 2 which is bounded away from zero and which is 
periodic in both x and y. For simplicity, we will assume that the basic cell of 
periodicity is [0, l] 2 , i.e. that £0(2: + 1, v) = £o( x , y) — Eo(x,y + l) for all (x, y) £ K 2 . 
For any open set O C K 2 we denote by L 2 g (0) the weighted £ 2 -space with norm 
given by 

IMleo = / £ o(x)\u(x)\ 2 dx. 
Jo 
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A 2D-periodic crystalline structure may be modeled by an operator Lq in L 2 Q (R 2 ) 
given formally by 

L u= t — rAu with domain D(L ) = H 2 (R 2 ). 

This is a self-adjoint operator. A standard tool to analyse periodic problems is the 
Floquet-Bloch transform. We will state here some of the results used in this paper 
and refer the reader to [TTJ [18] for proofs and more background on the theory. The 
Floquet-Bloch transform U x associated with the periodicity in the cc-direction is 

(1) 

U x : L 2 £o (R 2 ) -> L\ Q (n x (-7T, tt)), U x f(x, y, k) = ^y 1 ' 2 £ e* km f(x - m, y) 

where ft := (0,1) x R (more precisely, U x is defined by the formula (pj first for 
functions / with compact support, then extended to all of L 2 (R 2 )). 

We now consider an operator family L$(k x ) on the strip parametrised by 



k x G [— 7r,7r]. L a (k x ) is the self-adjoint operator in L 2 (Q) given by 



(2) L (k x )u= t^Au 

defined on the space of all functions u <E H 2 (£l) which satisfy the quasi-periodic 
boundary conditions 

(3) tt(l,i/) = e*-u(0,y) and y) = e ik ' g(0, y). 

It follows from the general theory [TTJ [TB] that Lq is the direct integral of the 
operators Lo(k x ): 



(4) L = J L Q (k x )dk x . 

[ — 77, 7r] 

As a consequence, the spectrum of the problem in the plane is 



(5) <r(L )= IJ CT ( L o(M). 

kx £ [ — 7T,7r] 

Moreover, in view of periodicity in the y-direction, similar arguments apply for 
each operator L (k x ) and the spectrum of the operator Lo(k x ) is 



a(L (k x ))= (J a(L (k x ,k)) 

fce[-7r,7r] 

where Lo(k x , k) is the operator — eo ^ x y ^ A on the unit cell (0, l) 2 with quasi-periodic 
boundary conditions in both the x- and y-directions: 

f u(l, y ) = e i *-ti(0 1 y) 1 f (!,») = e*-^(0,y), 



(6) 



u(ar, 1) = e* k u(x, 0), f^z, 1) = e lfe ^(x, 0). 



For notational simplicity we shall refer to the additional parameter as k rather than 
k y . 

From (J5j) , it is clear that any gap in the spectrum of Lq must arise from gaps in 
the spectrum of all the operators Lg(k x ). In this paper, we will fix k x and assume 
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that there are gaps in the spectrum of Lo(k x ), an arbitrary one of which we denote 
by (/io,Mi)- We will assume throughout that /io > 0. 

We now turn to the problem which is our main interest in this paper. We perturb 
the original periodic problem in the plane to consider a wave-guide W = R x (0, 1) 
inside the 2D-periodic crystalline structure. This new problem is modeled by an 
operator L acting on L^(M. 2 ) given by 

Lu = -Au 

e(x,y) 

with e{x,y) = Eo(x,y) + £\{x,y). The perturbation ej is supported in W , in the 
x-direction it is periodic with period 1 and it is such that inf eg + £i > 0. 

In view of the periodicity in the x-direction we can use the Floquet-Bloch trans- 
form ([I]) to generate, as in the case of L , a selfadjoint operator family L(k x ) on 
the strip fi, acting in the space L^(f2), parametrized by k x € [— 7r, it] and given by 

L(k, r )u := Au 

subject to the quasi-periodic boundary conditions ([3]). As before, the spectrum of 
the waveguide problem L is 



(7) <t{L)= |J a(L(k x )). 

k x £ [ — 7T,7r] 

The chief goal of this paper is to compare the spectra of Lq and L. In view of ([5]) and 
([7]), this amounts to comparing the spectra of Lo(k x ) and L(k x ). In [1] Lemma 10] 
it is shown that, a(L(k x )) can differ from a(Lo(k x )) only through the introduction 
of extra eigenvalues. Moreover, the essential spectra of the two operators coincide. 
In particular, the eigenvalues of L(k x ) cannot accumulate at any point inside the 
spectral gaps of L (k x ). 

We shall therefore study the eigenvalues of the perturbed problem 

(8) — Au — A(e + £i)u on ft 

where A € (/io,/ii), i.e. A lies in a gap of the spectrum of the operator Lo(k x ). It 
is understood that from now on all functions satisfy the quasi-periodic boundary 
conditions ((3]). From the equation ((8|) we get 

1 £\ 

Au — Au = A — u. 

It then follows that A is an eigenvalue in the gap iff 



(9) u = A (L (k x ) - A) 



-l / £i 



u 



holds for some non-zero u. Our strategy will be based on finding solutions of ([9]) 
using information on the resolvent of the unperturbed operator Lo(k x ). 

As mentioned before, we shall show that small perturbations of the operator 
Lo(k x ) create extra spectrum in an arbitrary fixed gap. Our second aim in this 
paper is to prove that the additional eigenvalues do not accumulate anywhere on 
the closure of the spectral gap (/kojMi), hi particular not at the endpoints fj,o an d 
Mi- 

We consider the non-accumulation result as a first step towards a better under- 
standing of the structure of the guided mode spectrum. One could surmise that 
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the guided mode spectrum a(L(k x ))\ a(Lo(k x )) can be written in terms of continu- 
ous band functions {Oj(k x )} depending on k x € [— 7r,7r]. However, since eigenvalues 
may be emitted and absorbed by the essential spectrum of the unperturbed operator 
Lo(k x ) as k x varies, the functions Oj may possibly be defined only on subintervals 
of [— 7T,7r]. An open question is whether a finite total number of band functions 
9j (possibly defined on subintervals of [— 7r, tt]) is sufficient to describe the guided 
mode spectrum. 

3. Eigenvalues at band edges and their eigenfunctions 

Our analysis of the waveguide problem will be via a study of the resolvent of 
the unperturbed operator Lg(k x ). This will be performed using the Floquet-Bloch 
transform. In this section we introduce a representation of the resolvent in terms of 
so-called Bloch functions and show some results which will be crucial in our later 
analysis. We note that the perturbation plays no role in this section. 

For fixed k x we consider the operator Lo(fc s , k) on (0, l) 2 introduced in the previ- 
ous section. The operator depends on k via the quasi-periodic boundary conditions 
and thus has a fc-dependent domain. We can transform the eigenvalue problem 
for Lo(k x ,k) into an eigenvalue problem for a fc-dependent operator with periodic 

boundary conditions in the y-direction: Let Vfc = V+i ( ^ ^ and define := V 2 , 
on (0, l) 2 subject to the boundary conditions 

f u(l,y) = e ik *u(0,y), f (1, y) = e ik * (0, y), 
{ u(x,l)=u(x,0), f!OM) = H(x,0). 
Then the self-adjoint analytic operator family — ^-A^ in L 2 o ((0, l) 2 ) is of type 
(A) (cf. [10]). i.e. the domain of the operators does not vary when k varies. 
Further, it is shown in |T0l Theorem VII. 3. 9] that there exist collections of functions 
{A s (fc)} se N and {</> s (fc)}seN which are real-analytic functions in the variable k on 
[— 7r,7r]. Moreover, for each s £ N, A s and </> s can be continued analytically to an 
open set 

{z e C : Re z e (-tt - S s , n + S s ), | Im z\ < r) s } 

containing the interval [— tt, tt]. For every k E [— tt,tt], the function (f> s (k) e 
L 2 o ((0, l) 2 ) is a normalized eigenfunction of — ^-A^ satisfying (|10l) and X s (k) is 
the corresponding eigenvalue. We note that the eigenvalues are not necessarily 
ordered by magnitude. Moreover, the analyticity results depend critically on the 
fact that we only let the scalar parameter k vary. The normalized eigenfunctions 
{ip a (x, y, fc)} se N of L (k x ,k) are then given by ip a (x,y,k) = e lky (f> s (x,y,k), with 
eigenvalues X s (k). We call these the Bloch functions. 

Lemma 1. Let \\u\\ 2 = ( ^ 2 \u\ 2 ) denote the unweighted L 2 -norm of a func- 
tion over (0, l) 2 . We have the following gradient estimates: 



(10) 



(11) ||V0 a || a < C7(v/A s (fc) + 1) and ||V^|| 2 < y/X 3 {k). 

Proof. As — j^Aj;0 s = X s (k)4> s (k), testing with cf> s , gives 

/ |V fc ^ s | 2 = / ( — -A k <j) s ) c/) s so = / A s (fc)0 s (/> s £ O = X s (k). 
J{o,i) 2 J(_o,i) 2 V £ o J J(o,i) 2 
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Hence, 

< ||Vjk0.|| 2 + \k\ Us\\ 2 < C(y/\ a (k) + l). 

2 

Note that, as k runs through a bounded set, the constant C can be chosen inde- 
pendent of k, but that it depends on Hl/eollno- The second statement follows from 
integration by parts: 

/ |V^.| a = / [-— A^ s )fe= / A s (fc)Vv0l£o = A s (fc). 
J(o,i) 2 ->(o,i) 2 V £ o ) J(o,i) 2 

□ 

We remind the reader that [1% is the lowest point of a spectral band and lies at 
the top end of a gap. The next result shows that only finitely many spectral bands 
can touch fi\ and that each A s (fc) can touch fii only finitely many times. 

Proposition 1. There are only finitely many pairs (sj,fej) such that X Si (ki) = fii. 

Proof. Let k, k G [— w,ir]. Arguing as in [TU1 VII. 3. 6], we get that there exist 
s-independent constants C\ , C2 such that 

\X s (k) - X s (k )\ < (Ci + |A s (fc )|)(e C2 l^l - 1). 

Next, choose \k — fcol < S, where S is such that e 2 — 1 < 1/2. Then 

\X s (k) - X s (k )\ < i(d + |A s (fc )|) 

and so whenever X s (ko) > we get 

X s (k) > ±X s (k )-C. 

For fixed ko, we have X s (ka) — > 00 as s — > 00. Hence 

inf X s (k) > -X s (kn) — C — > 00 as s — > 00. 

{fc:|fc-fc |<<5} 2 

Since 5 can be chosen independently of s and fc , we can cover [—it, tt] with finitely 
many intervals of length S and so mffc g r_ 7r)7r i A s (fc) — > 00 as s — > 00. Therefore, there 
are only finitely many values of s such that X s (k) — fix for some k. 

On the other hand, if for any fixed s we have X s (kj) = fi\ for infinitely many 
kj, then these must accumulate in [— 7r,7r] and the analyticity of A s would imply 
A s (fc) = fi\. But applying the well-known Thomas argument [19l [11], we see that 
none of the functions A s can be constant. □ 

We now discuss the expansion of functions in L 2 o ((0, l) 2 ) in terms of the Bloch 
functions ip s . Since for each k, the eigenfunctions ip s (', k) form a complete orthonor- 
mal system in L? E ((0, l) 2 ), we have for r G £ 2 (0, l) 2 and k G (— 7r,7r) 

r = Y^ rim, (r, if) s {-, k)) L 2 gmm t/>(-, k) 
and hence by using Parseval's identity and integrating over k, we get 

(12) IklLj ((o,i) a ) = 2^ / I ( r '^(-' fc ))L 2 o (( ,i) 2 ) I 2 dk - 

sen ^ 



\V<t>s 
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We now derive a representation for (Lo(k x ) — A)~V with r <G L 2 ((0, l) 2 ) in 
terms of Bloch functions. Whenever we apply an operator with domain L 2 (fl) to 
a function r £ L 2 o ((0, l) 2 ), we extend r to all of fl by zero. We will use the same 
letter for the extended function. 

Let U denote the Floquet transform in the y- variable defined analogously to (p} 
and set 

(13) P s (k,r) := <C/r(.,fc)> s (.J))^ ((01)2) ^ s (.,fc). 
We note that since r is supported in [0, l] 2 , we have 

(14) P s (k,r) = -^=(r,M;k)) L * ((01) M;k). 

V Z7T e u ' ' ; 

In view of the analyticity of the Bloch functions ip s , it follows that also P s depends 
analytically on fc in a small complex neighborhood of [— 7T, n] (note that k M> ip s (-, k) 
is analytic in fc, which can be seen by expansion into power series). 

Using this notation we have the following resolvent formula [TTJ [16j HB] f° r the 
operator L (k x ): 

(15) (L (fc x )-A)- X r= ^LV; / (A s (fc) — Xy 1 P s (k,r)dk 

for A outside the spectrum of Lo(k x ) and r g i 2 ((0, l) 2 )- 

4. Generation of spectrum in the gap 

In this section we use the representation of the resolvent by Bloch functions and 
the variational principle to show that if the perturbation e% is of fixed sign it will 
lead to the generation of extra spectrum in the gap (/xq,Mi) °f the spectrum of 
the operator Lo(k x ). For definiteness, in this section we assume that e\ is a non- 
negative function which is positive on a set of positive measure. Then there exist 
a > and a set D of positive measure such that infu e±\n = a. 

In set 

v — f^-u 
V £ o 

Then v is supported in [0, l] 2 , as £i|q is, and v satisfies 

(16) v = \[^- (L (k x ) - A)" 1 A 

V eo V £ o 

Note that vice versa, if v satisfies P^)l , then 

u := A (L (k x ) - Xy 1 J— v 
V £ o 

satisfies (|9|) and lies in L 2 o (f2). It is therefore sufficient for our purposes to study 
(HFJ. 

We now define the operator A\ on L 2 ((0, l) 2 ) by 



A x v = (xM(L (k x ) - Xy 1 x [^v) 

V V £ o V £ o J 



(0,1) 2 

and note that (|16| has a non-trivial solution if and only if 1 is an eigenvalue of the 
operator A\. 
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Lemma 2. Let A € 

compact. 



Then A\ : L\ ((0, l) 2 ) -» L 2 ((0, l) 2 ) is symmetric and 



Proof. Let w,i> € L 2 ((0, l) 2 ). Then 



(e A x u,v) L2aoi)2) 



eo* / — it, A 
1 V £ o 



-A - A 



£? 



L 2 (fi) 



L 2 (n) 



£ou, A* 



A- A 



(0,1) 2 / £2((0,1)2) 



= (£QW ) A«>i2((Q J l)2), 

so the operator is symmetric. Moreover, by standard estimates (see e.g. [3]) for the 
elliptic operator on the strip, 



-Ia-a 

£o 









< 


(-±A-*r\ 











ffi((0,l) 2 ) 

< CA ll U llL2 (n) = C\ ||w|| i 2(( 01 )2) ■ 

Thus A\ is the composition of a compact map with the continuous map of multipli- 
cation by the function and is therefore compact as a map from L 2 ((0, l) 2 ) — > 

L 2 ((0, l) 2 ). Multiplication by the bounded and boundedly invertible weight £o does 
not change this. □ 

We now investigate the dependence of the maximum eigenvalue of A\ on A. First 
define 

(A\u, u)c 

l \ \ \ ^ ' ' en 
Kmax(A) := SUp — r . 

l|„H/o {u,u) So 

Then if K rnax (X) > 0, it is the maximum eigenvalue of A\. 

Lemma 3. (1) On the interval (/^o,/-ti) the map A i-?> Kmax(A) is continuous 
and increasing. 

(2) If there exists A' € (/Lto,/zi) such that K max (A') > 0, then A H> K max (A) is 
strictly increasing on (A',/ii). 

Proof. (1) = A, /— (Lo(k x ) — A) -1 ./— , so A i-» is norm continuous as a 



map from (/io>Mi) to C{L EQ ) and for any A € (/io,A*i) an d an y £ > there exists 
5 > such that for |A — A| < 6 

(A^u,u) - (A A u,u) £ 



< Ua- A, 



lull 2 < £ ||u|| 2 
I 1 1 En — II lleo 



This implies 



K m ax{X) - K max (X)\ < £, so A n- K max (A) is continuous. 
Let /iQ < A < A < /ii. Then 

A A 



(A - A)A s (fc) 



A s (fc)-A A s (fc)-A (A s (fc)-A)(A s (fc)-A) 



> 
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since (X s (k) — X)(X s (k) — X) > and X s (k) > for all s and all k. Thus A H> K max (X) 
is monotonically increasing. This proves the first statement of the lemma. 

(2) We now investigate strict monotonicity of A H K max (A). We note that 
A i— > K max (A) need not be differentiable, as it is possible for eigenvalue branches to 
cross. However, K max is piecewise analytic on the interval [X',fii). Note that for 
A > A', K max (A) > K max (A') > by the monotonicity shown in the first part of 
the proof, so the point in the spectrum of A\ does not cause difficulties for the 
analytic continuation of the eigenvalues (the reader may refer to the discussion in 
Kato's book, [TUl Theorem VII.3.9 and Remark VII.3.11]). 

Consider now some subinterval of [A',/^i) where K max is analytic and let u\ be 
the corresponding normalised eigenfunction of Ax, depending analytically on A. 
Now, 

K max (A) = {Axu Xl ux) £g , 
so using the symmetry of A\ we have 



<9K max (A) 
dX 



dA x 
dX 

dA_x 
dX 



dA> 

dX 



, du x \ , I , du x 

Ax ^x> ux / e + \ AxUx 'Tx, e 



"A, 



dux \ 
dX A 



X s (k) 



A^(k)-xy 



— Ux,1pi 

£o 



dfc, 



where in the last step we have used the representation of the resolvent via Bloch 
functions (|T5j) . Now, X s (k) > for k ^ (else the operator L (k x , k) would have a 



non-positive eigenvalue). Assume 



OX 



<(A) 



0. Then 



for a.e. k 



and all s 6 N. Thus . —ux — which implies Axux = 0. However, this cannot be 



the (A) > 0. Therefore, K max (A) is strictly increasing on (A',^i). □ 



We are now ready to prove our main result which shows that small perturbations 
of eq create eigenvalues in the spectral gap. 

Theorem 1. Assume that 

(pi - jU ) inf e 



(17) lki|| 00 < 



Then there exists an eigenvalue of the operator L(k x ) in the spectral gap (/io,/ii) of 
L {k x )- 
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Proof. We shall use the Rayleigh quotient to obtain information on the eigenvalues 
of A\. Let u £ L 2 ((0, l) 2 ), then using (fTo)l we have 



(£ A x u,u) L2aoA)2) 



A ( e 



A- A 



£ 



J (0,1)2 



A( e [— A- A 



£l £i 

— u, a — u 

£o V £ o , 



L 2 (0) 



(18) 



J - 7T sen 



-u,ip s (-,k) 



dk. 



To first get an upper estimate on the Rayleigh quotient, let s' be such that /xi 
is the lowest point of the s'-band and no is the highest point of the (s' — l)-band. 
We note that such an s' must exist for there to be a gap. Let A £ (/zo,/xi). Then 
since X s (k) — A < for s < s' we have 



(e Axu,u) L 2 {(() 4)2) < — J ^2(X s {k) - A)" 1 

E 



— u,tp s (-,k) 



dk 



< 



< 



(19) 



2rr&*i - A) 
A 

2tt(aii - A) 
A 

(Mi - A) 



s>s' 



E 

sgN 



— U,1p s (-,k) 



-u,i> s {-,k) 



dk 



dk 



< A INL I, i|2 

" (/ii - A) inf £ " £0 



Therefore, if the perturbation ex is sufficiently small such that (fT?]) holds, we 
can find A' € (//o,Mi) such that 



(20) 



Kmax(A') = SUp 
II" II #0 



(A A 'U,m) 



< 1. 



We next want to find a lower bound on the Rayleigh quotient. Let kg be such 
that A s ' (feo) = Mi- 

We remind the reader that e\ > a > on a set .D of positive measure. It follows 
from unique continuation (see e.g. [17]) that for any Bloch function tp s (-,k) we have 



|^(-,fc)| 2 eo >0. 



Therefore, and by continuity of the Bloch functions there exist 6 > and a > 
such that 



I (ips>(-,k ),ip s >(-,k)) L 2 {D) | 2 > a 
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for fc e (fco — 5, ko + S). Now choose u = y fjV's'O, &o) 011 D and extend u by zero 
to (0, l) 2 . Then from JTSJ we get that for A e (Mo,Mi) 



A 



> 



> 



> 



2tt 


|| U 


2 




£o 




A 




2tt 


|| U 


2 




£o 




oA 




2tt 


\\u 


2 






£o 




1 




2tt 


|| U 


2 




£0 




a A 




2tt 


|| U 


2 




£0 



^(A^-A)- 1 (^(^oj.^t'^))^ (Z3) 



E 

77 s<s' 



(A.(fc)-A)- 1 (M;ko),M;k))w id) 



dk. 



dk. 



dk 



-5 K'(k)-X 
A 



^,A-A s( fc) 

ko+S 



dk 



dk 



— C\, 



k -s X s '{k)-\ 

where C\ is bounded as long as A stays away from fiQ. Using the Taylor expansion 
of the analytic function A S ', which has a minimum at fco, we see that 

Xs'(k) < Ml +a n {k-k ) 2n 

for fc e (fco — 5, kg + S) and for some n £ N with 



1 



Thus for 6 < 1 

1,1:0+6 dk 



fe -5 Ay(fc)-A 



> 



(2n)! {G[fe -.5,feo+'5] 



rifc 



A^(0>0. 



fco-5 Ml - A + a n(k - fco) 2 ™ 
1 



> 



k +S 



dk 



,«i - A 



■ arctan 



and hence 



(£0^AW,w) 



/ii - A 
arctan 



_ 5 /ii - A + a„(fc - fc ) 2 
/ii - A 



o(A)> 



/xi - A 
oo as A 



X = -5 

oo as A /ii 



Since by (|20[) we have some A' with n m ax(X') < 1 and K maK (-) is continuous, we 
can find A £ (A',/ii) such that K max {\) — 1. From ([§]) and (fll))) we see that this 
gives an eigenvalue of the perturbed strip operator L{k x ) between A' and Ml- □ 

Remark 1. (1) Our result only holds for small perturbations. One would of 
course generally expect that larger perturbations lead to the creation of more 
spectrum in the gap as lower eigenvalues of A\ cross 1, but our method of 
proof does not cover those cases. However, existence of guided modes for 
sufficiently large perturbations is shown for example in [51 [TS] . 
(2) We consider a special case: Assume that £i(x) = ai(x) for some function 
e(x) > and a £ K. + . This enables us to study the spectrum in terms of the 
scalar parameter a. We can interpret higher values of a as switching on or 
increasing the perturbation. Let A\ = A\/a and f ma x(A) be the maximum 
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eigenvalue of A\. Then we have A € a{L(k x )) iff cw m ax(A) = 1. This gives 
the following results: 

(a) Monotonicity ^max(A) and the fact that the maximal eigenvalue 
^max(A) of A\ tends to +00 as A — > fix implies that for any given small 
a, we can find A such that ^max(A) = — , so an additional spectral point 
is produced as soon as the perturbation is switched on. 

(b) Strict monotonicity o/A 4 ^max(A) implies that there exist A' > (mj 
and /3 > such that v max : (A',/11) — > (/3, 00) is invertible. Given 
a < 4, i/ie function g(a) := i/ ml ^(i) always gives an eigenvalue and 
g(a) — > jtti iti/ien a — > 0. 

(c) Let A G (/io,Mi)- TTie estimate (|19l) shows that the maximum eigen- 
value v max (X) of Ax satisfies 

v (X)< A " £ ~"~ 

- A)mf e 

Hence, if a < ^xjj^p^ 5 a^ max (A) < 1, so A is not an eigen- 

value. This implies that the perturbation needs to have a certain size 
before the spectrum can appear at any given point in the gap away from 

Ml- 

(3) In the case when E\ is a non-positive function analogous results can be 
shown by considering the operator 



A' x v = —xJ — — (Lo(kx) - A) 1 J— —v. 
V £ a V £ o 



In this case, the extra eigenvalues created for small perturbations will appear 
at the bottom end of the spectral gap near iiq . 

5. Non-accumulation of the eigenvalues 

In this section we wish to show that new eigenvalues generated by the perturba- 
tion cannot accumulate at the band edges. Unlike in section 21 we no longer make 
any assumptions on the sign of the perturbation. 

In ©, set 

v = — u. 

£0 

Note that v is supported in [0, l] 2 , as e%\n is, and v satisfies 

v. 



(21) v = \ £ -± (L (k x ) - Xy 1 : 



£0 

Therefore, we shall study the spectrum of the operator A— (Lo(k x ) — A) - acting 



on functions supported in [0, 1 



2 



5.1. Analytic continuation of the resolvent. Recall the resolvent formula (|T5|) : 
(22) ( Lo ( kx )-X)- 1 r = -^=y] f (A.(fc)-A)- 1 P,(fe,r)dfc 

for A outside the spectrum of L^(k x ) and r G L 2 (Q,) with compact support. Re- 
stricted to functions with compact support, the resolvent (Lo(k x ) — A)^ 1 is compact 
and we would like to use meromorphic Fredholm theory (see Theorem [2|) to analyze 
the spectrum. However, the resolvent is not well-defined in a neighborhood of the 
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band edge. To overcome this difficulty, we transform the critical integrals in (|22|) 
by integration over a suitable contour in the complex plane (see Figure [T]) instead 
of the real interval [— 7T, vr]. In this way, we obtain an analytic operator family, see 
(|23)) . that coincides with the resolvent on a sector of the complex plane. We will 
do this construction for the lower edge of a band near fi\. A similar construction 
is possible near the top end of a band. 

Note that by Proposition[T]there are only finitely many pairs (s p , k p ), p = 1, N 
such that for A = /Zi the integrand (X Sp (k) — A) -1 P Sp (k,r) in (1^1) is singular at 
k p , since X Sp (k p ) = n\. As none of the analytic functions A s can have a zero of 
infinite multiplicity by the Thomas argument (see [1HI19] ). we have for eachp that 
X^J(k p ) — for j — 1, ...,m p — 1 and xi" lp \k p ) ^ for some even m p > 2. Since 
Hi lies at the bottom end of the band, we can write X Sp (k) = fi% + (k— k p ) mp g p (k) 
with g p an analytic function of k, g p (k p ) > and <? P (K) C M. 

We set h Sp k p {k) :— (k—k p ) m %/ g p {k). Here, we choose the branch cut of the root 
away from the positive real axis, so that it does not intersect with the set {g p {k)} for 
k in a small neighborhood of k p . Since h' Sp k p (k p ) = mp / g P (k p ) ^ 0, we can locally 
invert /i Sp .fc p to obtain k = hj 1 k (v), and k will depend analytically on v for v in 
a ball M centered at 0. Let Af+=Nn {Im(v) > 0} and N~ = Af (1 {Im(^) < 0}. 
Clearly, h Spikp : R — > R, so h~* kp : A/' ± -)• C\R. As A/' ± are connected, K^i^) 
are connected and are therefore completely contained in either the lower or upper 
half-plane. Take v = ±is e A/' ± . Then, as 

^ fc J'(o) = 1 / m v / ^(M> ' 

by Taylor series expansion we get 



h s v \J ±is ) = k P+ I 7=^F^ + r ±( is ) 



±is 
V5p( fc p) 

with r±{is)/s — > as s — > 0. Hence, lni(hj 1 k (is)) > and lm(hj 1 k {—is)) < 0, 
so hj 1 ^, maps J\f + to the upper and N~ to the lower half-plane. 

Proposition 2. All solutions of (k — k p ) mp g p (k) — v mp in a neighborhood of k p 
are given by 

k = K\ (e^/^v) 

Sp , Kp \ ' 

with p = 0, m p — 1. 

Proof. In a punctured neighborhood of k p we have the following equivalences: 

«*■ . =e - 2 "P/m P) p£ {Q,...,m p - 1} 

(fc - k p ) m K/g P (k) 



fc = ^\ (e 2 " p/m ^), p G {0, m p - 1}. 



□ 



For simplicity of notation, from now on we restrict ourselves to the case when 
only one band, which we call the so-band, touches fix. All results generalise in the 
obvious way in the case when more than one band touches 
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Let (so,fc p ), p = 1,...,N be all pairs described in Proposition Q] with m p being 
the order of the first non- vanishing derivative of A So at k p . Let m be the lowest 
common multiple of the m p and q p — m/m p . Using the Taylor expansion of X So (k) 
around each of these k p , we find a complex neighborhood N(fi\) of \x\ and balls 
BrQci), Bn(kisi) of some radius R around each of the k p such that by Proposition 
[2] the equation A So (fc) = /ii + /i m " has precisely m p solutions in Bn(k p ) whenever 
fii + € N(fix). Moreover, we can find a smaller neighborhood N(/j,i) of /ii 
such that for each i = 1, N and all fii + /i mp € N(/j,i) all solutions in Bn(k p ) of 
A So (fc) = /ii + /i mp in fact lie in a ball B R / 3 (k p ) of radius i?/3 around 

We now choose a contour G as indicated in Figure [1] 




Figure 1. The contour G. 

For later estimates, it will be useful to have a bound on the distance from the 
curve G to solutions of the equation A So (fc) = z. This is given by the following 
lemma. 

Lemma 4. There exists another neighborhood N C N(fii) such that for some 
positive number Sq we have dist(G, {k € C : A SQ (k) — z}) > So for all z € N . 

Proof. For a contradiction, we assume that for all neighborhoods N and for all 
8q > there exists z G N such that dist(G, {k : A SD (fe) = z}) < S . This implies 
that there exists a sequence (z n ) in iV(/xi) such that z n — > fix and 

dist(G,{fc : A SD (fe) = z n }) -> 

as n — > oo. This in turn implies the existence of a sequence (fc„) such that A So (k n ) — 
z n and dist(G, {k n }) ~ > as n — > oo. 

Assume for another contradiction that fc„ e Bn(k p ) for some i. Then by defi- 
nition of N(fxi) we would have fc„ G B R / 3 (k p ) for large n. However, this contra- 
dicts that the k n are very close to the curve G for large n. We can deduce that 
k n Bfi(k p ) and since dist(G, {k n }) ->0we get that Im k n —> as n —> oo. 
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Now, [— 7r, 7r]\ U (Bn(kp) n R) is compact, so there exists a subsequence k n . 
converging to some k £ [— 7T,7r]\ U (-Bft(fcp) n R). It then follows that A So (fc n ) — ¥ 
X So (k), and z n . Ml- So X So (k) — Mi- Hence, k is an additional real solution to 
A So (k) = /ii contradicting the assumption that the k p are all solutions. □ 

We are now in a position to introduce the operator family that will provide an 
analytic extension of the resolvent. Let O be a small neighborhood of such that 
Mi + /i m € N for any fi e O. For \i e O \ {0} and r E L 2 eo {{0, l) 2 ) we define the 
following operator family 



B^r) = 4=^ r(A.(fc)-A*i-/* m r 1 i'.(fc,r)d* 

V27T , J-7T 



1 



(23) + -L / (A S0 (fc) - mi - M™) -1 A (*, r)d* 



G 
V 

2^E £ ^ o (^ fci (e 2 " p / m4 /i 9 0)- 1 ^a(^^Je 2ffWmi M 94 ) ) r-). 

i=l p£Pi 

where, P,; = {0, ^ — 1}. The next result shows that does indeed extend the 
resolvent. 

Proposition 3. For \i 6 \ {0} with axgfi € (0, 2-K/m), we have 

B^r) = (L {k x )-^ 1 -(i m )- 1 r 

/ O rr€L 2 Q ((0,l) 2 ) 

Proof. Noting the resolvent formula (fl5t . we see that £? M is obtained from the 
resolvent by replacing integration over (— n, n) with integration over G when s — 
So and adding the final term with the point evaluation. Hence, this amounts to 
a simple application of the residue theorem: We need to find the solutions to 
X So (k) = Mi + M TO ; i- c - to (k — ki) mi g(k) = fi m , lying between G and the real axis. 
Since h~ x k . maps points in the upper half-plane to points in the upper half-plane 
and points in the lower half-plane to points in the lower half-plane, by Proposition 
[2] and the discussion following it, these are given by k = hj^ k (e 27T ' lp ^ mi /j, m ^ mi ) for 
those p € {0, rrii — 1} such that e 27n P/ m ifj l m / m i \[ es m the- upper half plane. Now 
for arg/x £ (0, 2-K/m), 

arg ( e 2^p/m ifJ rn/mA g ( 2 ^V_^V_ + . 

V / \ rrii rrii m, J 

This is in (0, tt) if and only if p € Pi. Note that in calculating the residue, we have 
used that Res z f/g = f{z)/g'(z) whenever g(z) = 0, g'(z) ^ 0. □ 

5.2. Properties of the operator family B^. We next look at properties of the 
operator family B^ from ((23|) with m in the small neighborhood O of 0. Recall that 
our aim is to be able to apply Fredholm theory to this operator family. 

Proposition 4. fi n- B^ is analytic as a map from O \ {0} to £(L 2 o ((0, l) 2 )). 

Proof. We discuss the three terms on the right hand side of (J23J) separately. 
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(1) The only problem in the residue term arises when X' Sg (h s ^ k . ( e 27 ™p/ m i^<3; )) 
0. Now, /i = is locally the only solution of h~ 1 k . (/i) = fcj and, since 

A' S0 (fc) = (fc-fc i r'- 1 [m ig (fc) + (fc- W(fc)] ; 

for fe=fei 

fcj is locally the only zero of k H> A' s (fc). Therefore, the residue term is analytic for 
veO\{0}. 

(2) We next consider the contour integral 

(24) -±= f (X ao (k)-m-ii m )- 1 P^(k,r)dk. 



/2n 

Lemma 2] implies that along the contour G we have that |(A g (fc) — X so — /i m ) _1 | is 
bounded. In particular, this implies that the integrand in (IM1) is analytic in /i. 

(3) Finally, we consider the infinite sum component in the resolvent. We note 
that by an argument as in the proof of Proposition!!] there exists rj > such that for 
all band functions X s which do not touch /ii and for all k, we have |A s (fc) — yLtx ] > V- 
So there exists a neighborhood O of the origin such that for all /z € O and for all non- 
touching bands s we have \X s (k)—fii—fi m \ > | for all k. Hence, \(X s (k) — fii— ii m )~ 1 \ 
is uniformly bounded in s and /x. We need to test with a function <fi and show that 
for any (j) the resulting function is analytic in \i. Fubini's Theorem allows us to 
interchange the integration in k and the L 2 -scalar products. Then we have for 
H E O and M > N sufficiently large, 

g L a, w -mi-^ {Ur{ -> k) > *»« ^ kU) ^ dk 

< E r -\(Ur(;k),M;k)) £0 \\(M;k),4>) £0 \dk 

s=N J - 7! V 

2 



< 



^ / |([/r(,fc),^(,fc))|2 o dfc £ / 

\ s=N J - n \ s=N J - 



(V>s(-, k),4>) \ 2 eo dk. 



By completeness of the Bloch functions (see (|12p ). this tends to as N,M — > oo. 
Hence we get uniform convergence of the series in fi and the limit function is 
analytic. □ 

Proposition 5. is compact for /i £ 0\ {0}. 

Proof. Again, we consider the three terms on the r.h.s. of (|23[) separately. 

(1) The term from the residue is clearly compact, as it is a finite rank operator. 

(2) We next consider the contour integral (|M|) . Lemma 2] implies that along the 
contour G we have that |(A S0 (A;) — /Ui — (i m ) \ is bounded and the same is true of 
IIV'soO' ^)ll- Using (|14[) and taking the norm into the integral, we can now estimate 

||/ G (A S0 (fc)- / i 1 - A1 ™)- 1 P ;50 (fc,r)(.)*|| Hl((0!l)2) 



< CJ G p^ggj^ l (r,^ S0 (-M S0 \dk <C\ 
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where the second inequality follows using (fTT|) . Hence, as an operator in L^ o ((0, l) 2 ), 
this part of the resolvent is compact. 

(3) Finally, we consider the infinite sum component in the resolvent. We note 
that, as before, there exists 77 > such that for all fi € O and for all non-touching 
bands s we have \\ s (k) — fit — fi m \ > \ for all k. 

Hence, |(A s (fe) — fit — fi m )~ 1 \ is uniformly bounded in s and fi and 



- 1 P s (fc,r)(:)& 



M „7r 

£ / [X s (k)-fit-f 

s=N J - yT 

/ Yl ( A «( fc ) - A*i - M" 1 )" 1 (EM-, k), ^(-, fc)) £o ^(:, fc)dfc 



: C 



: C 



• C 



s=N 
M 



E( A *( fc ) 



mi 



Ai " l )- 1 ((7r(.,fc),V s (-,fc)) £n ^(:,fc) 



s=JV 
,1/ 



ff 1 
2 



rffc 



i? 1 



E 

M 

E 

s=JV 



l^0,fc)IUi 



(A s (k) - fit - fi m f 

2 



(Ur(;k),M;k)) en \ 2dk 



+ 1 



(X s (k) - fit - fl m )' 



:\(Ur(;k)^ s (;k)) e \ 2 dk 



M „7r 

< CE / |(C/r(,fc)> s (-,fc)) £ J 2 dfc^O 



as M, JV — > 00. In the third line we have used the H 1 -orthogonality of the eigen- 
functions {ip a (*, k)} se ^ for fixed k. Thus, we have a Cauchy sequence in -ff 1 ((0, l) 2 ). 
Moreover, the i/^-norm of the limit is bounded by C ||r|| e which gives compactness 
as an operator in L 2 Q ((0, l) 2 ). □ 

Proposition 6. has a pole of finite rank for fi = 0. 

Proof. We note that the only pole comes from the residue term at fi = 0. The 
operator has a pole of finite order as A' so only has simple zeroes of order rm — 1 
at the h and hj 1 ),. is analytic in a neighborhood of 0. Moreover, the factors 

p so{KlkM 27Tip/m '' i 'v qi )> r ) are 01 rank L D 

Proposition 7. Le£ fi E \ {0} suc/i i/iai arg/x G (0, ir/m). Let 



(25) 



£0 



T7ien (/ — -B/j)^ = only has the trivial solution in L 2 Q ((0, l) 2 ). 

Proof. Suppose that there exists v ^ such that (I — B^)v = 0. Set w :— 
(Lf)(k x ) — A) X w where A = fit + /i m . Then the equation for v implies that 
v = X—w, so (— — A — A) w = v = X—w, or —Aw = A(e + £i)w. But this 



implies that A = fit + fi m is a non-real eigenvalue of the operator L(k x ), yielding a 
contradiction. □ 
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5.3. Main result. To prove our result on the non-accumulation of eigenvalues, we 
make use of the following theorem on meromorphic Fredholm theory |18[ Theorem 
XIII. 13]: 

Theorem 2. Let D be a domain in C, S a discrete subset of D. and H a Hilbert 
space. Assume we have a family of operators {A z : z £ D} such that 

(1) z t-¥ A z is analytic as a map from D \ S to C(H), 

(2) A z is compact for z £ D \ S . 

(3) A z has poles in S of finite rank, 

(4) there exists z 6 D\S such that (I — A z )u = has only the trivial solution. 

Then there exists at most a discrete set S C D such that (/ — A 2 ) _1 € C(H) for 
allzeD\(SUS). 

We can therefore prove our main result of this section: 

Theorem 3. The spectrum of the problem on the strip, i.e. the spectrum of the 
operator L(k x ) can not accumulate at the ends of bands of the spectrum of the 
operator Lo(k x )' 

Proof. By [TJ Lemma 10], the spectrum of the operator L{k x ) outside the bands 
can only consist of eigenvalues. By (|2"Tj) . it is clear that A = fj,\ + /i m <E a p (L(k x )) 
if and only if / — is not invertible. 

We apply Theorem [5] to the operator family from (|25l) with /i in a small 
neighborhood O of 0. The properties (l)-(3) follow from the propositions proved 
for B^ in the previous sub-section which obviously carry over to B^, while (4) is 
shown in Proposition [7] □ 

Acknowledgements. The authors would like to thank T. Dohnal, who showed 
them a numerical example illustrating the results of section 4. 
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